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(54) Method and apparatus for calibration of a multi-channel fiber optic interferometric sensor 
system In a signal processing system 



(57) A calibration tecfinique tor lasers of a data ac- 
quisition module in a signal processing system, wherein 
an optimal value of a phase modu lation coefficient is cal- 
culated for each laser and used as the amplitude of a 
cosine wave signal for the laser module carrier frequen- 
cy for phase modulation of the laser output signal. The 
phase nrKxiulation coefficient Is determined from a gb- 



bal average l/Q ratio, calculated for all channels of the 
selected laser module, as an average of all i/Q ratios of 
in-phase and quadrature phase angle components of 
the signals returned from the sensors. The phase mod- 
ulation coefficient for the selected laser module is de- 
termined empirically or by dividing the previously deter- 
mined phase modulation coefficient by the previously 
detennined global average i/Q ratio. 



F/Gf 



104 



CM 

< 

CM 
ID 

s 



Q. 

LU 



L! 



ARRAY 1 I ARRAY | | ARRAY | 
PANEL 1 I PANEL i I PANEL | 
















ARR 
MNf 


AY 

1 


ARRAY 
PANEL 


ARRAY 
RAMEL 



120 

/ 

X(#(t)) 



20^aOB^ 



INBOARD RECEIVER 



fie 



LASER MODEL 125 



PHASE MODULATOR 



U9 

LASER DRIVE CARD V^T 



LASER CONTROLLER 
CARD 



i46^ 

AMP/FREO 
CONTROL 

TIMING 



riie 



V" 

INBCARD _ 
RECEiVER POW^ 



'f/2 



POWER 
DISTRIBUTION/ 
CONDITIONING 



FREQUENCY 
SYNTHESIZER 



f40 



MEMORY 



/ Jg ^CM P/STATUS 



106 



939 



CPU 



CMD/ STATUS 



RECEIVER 
-^8D 



216- 



CONTROLLER^ y 



POD/|| ADC/llDigital 
RCVR|Digi1al bemod 



(DemuK Hsignol 



VME 
I/F 



Roceway 
I/F 



VME 
I/F 



>fS4 



t24 



f26 



NO 



COMMON 




FIBER 


VTAMFQRACR 




CHIV^fCL 


CARD 




CARD 



Printed by Jouve. 75001 PARtS (FR) 



BEST AVAILABLE COPY 



EP0 936 452 A2 



Description 

[0001 ] The present invention relates generally to the 
field of signal processing, and nriore specifically to a 
method and apparatus tor calibration of a multi-channel 
fiber optic interferometric sensor system in a signal 
processing system. 

[0002] Acoustic listening systems for undenwater ap- 
plications are well known in the art. They are typically 
equipped with arrays of acoustic sensors which provide 
sensitive undenwater listening capabilities and may pro- 
vide relative position information. Each sensor reacts to 
an incoming pressure wave by modulating an input sig- 
nal and the outputs of all the sensors are processed to 
determine the sound and position information. These 
sensor arrays are generally either attached to the sub- 
marine hull or are towed behind the submarine. Ideally, 
the sensors would be attached to the submarine hull, 
but prior art acoustic sensors were too heavy for many 
submarine applications. Recent advances in acoustic 
sensor array technology, however, have produced ar- 
rays which are light enough to be mounted on a subma- 
rine hull and which still provide very sensitive signal re- 
sponse. This weight reduction has also albwed the 
number of sensors to be increased. 
[0003] Another application for undenvater acoustic 
sensors is in the geological survey industry, specifically, 
for underwater oil exploration. Vast arrays of sensors 
may be placed on the ocean bottom in the vicinity of 
Icnown oil reserves. A surface ship initiates an acoustic 
pressure wave with a large air burst. The acoustic pres- 
sure wave and its reflection off the ocean floor are de- 
tected by the sensor arrays. The data from the sensors 
is then processed and analyzed to determine optimum 
drilling locations or to monitor the status of known res- 
ewoirs. 

[0004] As the number and complexity of the acoustic 
sensors in these and related applications have in- 
creased, the associated signal processing electronics 
modules have likewise increased. Typical conventional 
systems use analog circuitry to interrogate the sensors, 
but these analog systems are subject to drift and it is 
very difficult to accurately calibrate them. Recent devel- 
opments have provided fiber optic interferometric sen- 
sors, as described in the article "Homodyne Demodula- 
tion Scheme for Fiber Optic Sensors Using Phase Gen- 
erated Carrier" by Anthony Dandridge, Alan B. Tveten, 
and Thomas G. Giallorenzi, IEEE Journal of Quantum 
Electronics, Vol. QE-18. No. 10, October 1982, incorpo- 
rated herein by reference. This reference teaches a 
modulation/demodulation technique in which the fre- 
quencies of I and Q components are different, which is 
employed in the present invention. 
[0005] The above-mentioned related patent applica- 
tions describe techniques for interrogation of fiber optic 
interferometric sensors in multi-channel applications 
and for fine tuning of a receiver card 122 by calibratwn 
of digital down converters of the signal processing sys- 



2 

tern 100 of Figure 1 . The signal processing system 1 00 
is comprised of a source end with lasers, sensors and 
control electronics, and a receive end. The sensors are 
located In-between the source end and receive end as 
6 part of the outboard telemetry. The sensors are fiber op- 
tic sensors which get the information from the environ- 
ment sensed by the fiber optic interferometers. The op- 
tical fiber interferometer has a reference fiber and a sig- 
nal fiber. A relative optical phase shift between these 
10 two fibers is the information returned by the telemetry 
and extracted by the receive end. 
[0006] Figure 1 is a high-level block diagram of the 
signal processing system 1 00. An inboard receiver mod- 
ule 102 containing the requisite electronic control and 
15 processing circuitry is located inside the submarine. A 
lightweight array of fiber optic interferometric sensors 
1 04 is located on the outside hull of the submarine. Each 
sensor in the array reacts to an incoming acoustic pres- 
sure wave by phase modulating a light signal 118 sent 
20 from a laser module 11 6 of the inboard receiver module 
1 02. Modulated signals from different sensors, each op- 
erating at a different carrier frequency are passively 
multiplexed onto a single fiber optic cable return signal 
1 20 and sent back to the receiver 1 02. The channel sig- 
2S nals are then demultiplexed and demodulated by the re- 
ceiver 102. 

[0007] The inboard receiver 102 is connected to the 
submarine's power system (not shown) via a power con- 
nection 106. The inboard receiver 102 has a power dis- 

30 tributlon and conditioning block 1 12 to provide the nec- 
essary power requirements for the receiver compo- 
nents. Each input light signal 118 is generated by the 
corresponding laser module 116. Each laser module In- 
cludes a laser 117, a phase modulator 119, a laser driver 

35 card 97 and a laser controller card 98. The laser 117 
used in the preferred embodiment is a f^odel 1 25 200 
mW Nd:YAG laser, available from Lightwave Electronics 
of Palo Alto, California. The phase modulator 1 1 9 is pref- 
erably a lithium niobate phase modulator. 

40 [0008] As shown in Figure 2, eight laser modules with 
eight lasers 117a - 117h drive eight corresponding sen- 
sor arrays 104a - 104h. Each laser 117 powers 56 sen- 
sors (channels), typically 28 channels on each side of 
the submarine. The laser 117 crystal is heated with the 

45 temperature level provided from a frequency synthesiz- 
er 142 (Figure 1). The frequency synthesizer 142 con- 
trols the temperature of each laser 1 1 7 crystal separate- 
ly, as is well known in the art. to produce laser outputs 
with eight different "colors" i.e.. different wavelength 

so nominally at 1319 nanometers, which are provided to 
the corresponding phase rTK>dulator 119a - 11 9h. 
[0009] On their way back to the receive end of the in- 
board receiver 102. all first channel outputs from each 
sensor array 104 are multiplexed together via a signal 

55 multiplexer 121. This is repeated for alt the other 55 
channels. Therefore, 56 fiber optic cable return signals 
120 return from the sensor array 104, with each fiber 
carrying eight multiplexed signals. The multiplexed sig- 
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nal is then prcxessed by a receiver card 1 22 (Figure 1 ). 
Depending upon the number of sensors used, several 
receiver cards 122 nnay be needed. 
[001 0] The frequency synthesizer 142 contains the 
master system clock and controls the operation of the 
laser module 116 by defining laser parameters, such as 
laser color, power, and temperature, provided to the la- 
ser module 116 on an amp/freq control line 146. The 
timing information is provided to the receiver card 122 
via a signal timing line 144. A CPU 130 and an associ- 
ated memory 1 40 provide system level control and sta- 
tus information to the component of receiver 102 via da- 
ta linos 132. 138. 141. The CPU 130 has a VME inter- 
face 1 39, attached to a VME bus 1 38, which exists be- 
tween the CPU 1 30 and the receiver card 1 22. The CPU 
130 is connected to the frequency synthesizer 142 via 
a data line 141 

[0011] The receiver card 122 demultiplexes and de- 
modulates the retum signal 120 and outputs a signal to 
a beamformer card 1 24 which provides additional signal 
processing. A fiber channel card 1 26 provides the signal 
information to an external fiber bus 110. Each receiver 
card 1 22. processes the signals from 7 fiber return sig- 
nals 1 20, wherein each fiber has eight laser channels, 
or a total of 56 channels from eight different sensor ar- 
rays 104. 

[0012] Figure 3 is a detailed block diagram of the re- 
ceiver card 122 shown in Figure 1. The return signal X 
(0(t)) 1 20, containing 8 channels, output by the acoustic 
sensor array 1 04 is input into a polarization diversity de- 
tector (PDD) 200. An example of a PDD which has three 
outputs IS described in U S Patent No. 5,448,058, enti- 
tled "OPTICAL SIGNAL DETECTION APPARATUS 
AND METHOD FOR PREVENTING POLARIZATION 
SIGNAL FADING IN OPTICAL FIBER INTERFERO- 
METRIC SENSOR SYSTEMS." In the preferred embod- 
iment, a bi-cell or a two-output PDD is used. The PDD 
200 and some other elements of the receiver card 122 
prevent polarization signal fading in the return signal 
120. The PDD 200 converts the photonic energy of the 
optical fiber return signal 1 20 into two separate electrical 
currents 200a, 200b via two photodiodes. 
[0013] In each opto-receiver 202a. 202b a transim- 
pedance amplifier converts the input current to a volt- 
age, and a variable gain amplifier sets the voltage level 
in order to maximize the signal-to-noise level and to en- 
sure that the voltage level is below the saturation level 
of ADCs 204a. 204b. Anti-aliasing filters in the opto-re- 
ceivers 202a, 202b are used to filter the signal which 
then passes to the ADCs 204a, 204b, where the analog 
outputs of the opto-receivers 202a, 202b are digitized 
by high speed (> = 25.6 MSPS). high resolution (> = 
12-bit) analog-to-digitai converters (ADCs) 204a. 204b. 
The ADC used in the preferred embodiment has the part 
number 9042, and is manufactured by Analog Devices, 
Inc. 

[0014] The number of ADCs used is directly depend- 
ent upon the number of outputs from the PDD 200. Typ- 



ically, there are two outputs, but three are possible. At 
this point, the digitized output contains a complex signal 
comprised of all channels on the return optical fiber re- 
tum signal 120. which are the frequency division multi- 
s piexed phase generated carriers with the information- 
carrying sidebands. 

[0015] Each ADC 204a, 204b output is buffered and 
passed through to a multi-throw, muiti-pote bus switch 
206 which taps into each ADC 204a. 204b signal path. 

10 One output of the bus switch 206 goes to a digital down 
converter (DDC) 228, via signal path 238, which is em- 
ployed as part of a calibration channel, as described bo- 
low. The outputs of the bus switch 206 go to the signal 
channel DDCs 208a, 208b, 208n. The purpose of the 

IS bus switch 206 is to allow each of the DDCs 208 to be 
able to connect to any ADC output. This is required by 
the PDD 200 selection algorithm, discussed below with 
reference to Figure 7. 

[0016] Each DDC 208a. 208b. 208n acts as a digital 

20 demultiplexer by performing mix down and filtering of 
the digital information, separating one channel out of the 
composite signal. The DDC chip employed in the pres- 
ently preferred embodiment is the GC4014 chip manu- 
factured by Graychip, Inc. This particular device has 

25 2-channel capability. In other embodiments, each chip 
may only have 1/2 or 1 -channel capability. The number 
of DDCs needed depends upon the number of channels 
used in a given application. For example, if the return 
signal 1 20 has 8 multiplexed channels, then 4 2-channel 

30 DDC chips are needed. In the present embodiment. 28 
2-channel DDCs (56 channels) are used per receiver 
card 1 22, but only one calibration DDC 228 and receiver 
card controller 226 are required, as described betow. 
Thus, seven channel groups are needed to process all 

35 56 multiplexed signals, with 8 channels per multiplexed 
array output signal. 

[0017] Each DDC 208a, 208b, 208n outputs both in- 
phase (t) and quadrature phase (Q) data words that rep- 
resent the rectangular components of the phase signal 
40 from the bus switch 206. These I and Q components 
from the DDCs 208a, 208b, 208n are time division mul- 
tiplexed (TDM) onto separate I and Q buses. The I and 
Q components which define the same phase angle do 
not have the same frequency in the present invention. 
45 They are harmonics of the same frequency and one of 
these components is usually on the fundamental fre- 
quency and the other is on the first harmonic frequency. 
The output signal paths 210, 212 (which may be serial, 
unidirectional data paths or bi-directional data paths) 
so output the I and Q 8-bit data words to a buffer 214 which 
buffers the words and also converts the 8 bit words into 
16 bit words. 

[001 8] The I and Q words from all of the DDCs 208 in 
a time sequence are then demodulated. The I and Q 
55 words are converted from Cartesian to polar form by the 
coordinate transformer 21 6 to obtain the magnitude and 
the phase angle in polar coordinates. The coordinate 
transformer 216 of the preferred embodiment is a 
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TMC2330A chip manufactured by Raytheon Corp., but 
may be any similar device. 

[0019] The output of the coordinate transformer 216 
is the instantaneous phase angle, whose change is di- 
rectly related to the change in acoustic pressure of the 
acoustic signal from the environment, and whose rate 
of change Is directly related to the frequency of the 
acoustic signal from the environment for a given fiber 
optic sensor. The channel-by-channel Instantaneous 
phase angles are further processed in an acoustic signal 
processor 218, with integration and filtering algorithms, 
which is also attached to the VME bus 1 38. The output 
of the acoustic signal processor 218 may then be pro- 
vided to a visual display an oscilloscope, or to some 
other signal processing blocks, as desired. Two DACs 
222, 224 provide I and Q signals which can be used tor 
testing or viewing the output signal 
[0020] The calibratbn procedure described in the 
above-mentioned patent application U.S. Serial No. 
09/021 ,557 uses the receiver card 1 22 fine tuning loop 
for calibration of the receiver card 1 22 components. For 
that purpose the output 238 of the bus switch 206 is used 
to feed the calibration channel DDC 228. If there is more 
than one PDD 200. there will be more than one bus 
switch 206 feeding the calibration channel DDC 228, in 
a TDM manner. The calibration channel DDC 228 is set 
to mimic each signal channel in tum. The calibration 
channel provides a way to check PDD output selection, 
phase shift adjustment and i/Q balance on a per signal 
channel basis in a non-obtrusive way. For that purpose, 
the data is processed by the DDC 228 into I and Q words 
similar to the output of the signal channel DDGs 208a, 
208b, 208n. The I and Q words are then processed by 
the receiver card controller 226 (a digital signal proces- 
sor with its own RAM and ROM 232 memory) to extract 
the information necessary to determine channel per- 
formance. The receiver card controller 226 of the pre- 
ferred embodiment is a ADSP 2181 manufactured by 
Analog Devices, Inc. 

[0021] The operation of the receiver card controller 
226 is described with reference to the flowchart of Figure 
4. The following steps are performed by software which 
is either stored in the system controller's ROM 232 or is 
downloaded by the system CPU 1 30 to the signal proc- 
essor 226 DSP RAM upon system initialization. Upon 
power initialization or a system relnitializatbn 300, all 
interrupts are disabled at step 302. The receiver card 
controller 226 digital signal processing unit (DSP) is in- 
itialized at step 304, along with any associated RAM 
memory and the DDGs 208. The interrupts are then en- 
abled at step 306. A signal channel and a calibration 
channel are selected at step 308. Then, the best PDD 
signal is selected at step 310 by making the ADC selec- 
tion for each signal channel. This PDD selection step is 
further described below with reference to Figure 7. 
[0022] A phase shift algorithm for the in-phase I com- 
ponent is performed at step 31 2, which is shown in detail 
in Figure 5. Similarly, a phase shift algorithm for the Q 



component is performed at step 31 4. as shown in Figure 
6. An I/Q balance algorithm is executed at step 316, 
which is shown in detail in Figure 8. Steps 308 - 316 are 
repeated for each channel, until all the channels have 

5 been calibrated. In one embodiment, calibration Is per- 
formed only at system start-up. In the preferred embod- 
iment, however, the calibration procedure continues for 
each DDC for as long as the system Is operating. This 
results in a dynamic calibration system which unobtru- 

10 sively calibrates the system while signal processing 
functions are unaffected. 

[0023] The procedure for selecting the best PDD sig- 
nal (step 310) is illustrated by tho flowchart of Figure 7. 
At step 602, the first ADC 204a is selected, and I and Q 

15 data samples are obtained at step 604. A maximum 
'peak-io-peak value Is Calculated at step 606, when Q val- 
ue is near zero. When 0=0, the Ipeak-tc^peak "^^s maxi- 
mum value but it is difficult to catch that moment. There- 
fore, a set of I values is collected of points within an in- 

20 ten/al around the time when Q=0. There are preferably 
32 collected I values when the Ipeak-to-peak values are at 
their maximum Ipeak-to-peak value. These 32 maximum 
Ipeak-to-peak values are averaged to determine the aver- 
age maximum Ipeak-to-peak va'ue. Similar process may 

25 be performed to detemriine the average maximum 
Qpeak-to-peak value, When needed. 
[0024] The average maximum Ipeak-to-peak value de- 
termined at step 606 is then stored at step 608. Then 
the second ADC 204b is selected at step 610. I and Q 

30 data samples are then obtained for the second ADC 
204b at step 612. Again, an average maximum 
Ipeak-to-peak value is calculated at step 614. and this val- 
ue is stored at step 616. The two stored average maxi- 
mum Ipeak-to-peak values are then compared at step 61 8. 

35 The ADC 204 which produces the larger Ipeak-to-peak val- 
ue is then selected (steps 620, 622). Since selecting dif- 
ferent ADCs 204 will affect the I and Q signal levels sym- 
metrically, only one signal (I or Q) needs to be examined 
in this procedure. The PDD 200 selection step is per- 

40 formed for each channel separately. Thus, different 
channels within DDC 208 may actually be using different 
ADC 204 inputs, which is why each DDC 208 needs to 
be connected to both ADCs 204. 
[0025] The phase shift algorithms for the I component 

45 (step 31 2) and the Q component (step 31 4) are shown 
in Figures 5 and 6, respectively. The purpose of these 
algorithms is to compensate for the misalignment of the 
return signal 120 from the laser modules 116. and the 
internal receiver card 122 local oscillator signal in the 

50 DDC 208. The misalignment occurs due to the fact that 
the laser modules 116 and the receiver card 122 are all 
run by the same master clock from the frequency syn- 
thesizer 142. The adjustment is performed in the DDC 
208. The algorithm determines the phase offset values 

55 that result in a maximum phase signal (I data) and quad- 
rature phase signal (Q data) for each channel. This is 
accomplished by successively programming the cali- 
bration DDC's phase offset control registers with equal 
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phase increments and reading the corresponding I and 
Q data from the channel output registers. The i and Q 
data generated by the calibration DDC 228 are read by 
the receiver card controller 226 and stored into two sep- 
arate buffers. Once a sufficient number of data samples s 
has been accumulated, the phase shift algorithm is per- 
formed on the I and Q data. The outputs of the phase 
shift algorithms are an l-phase offset and a Q-phase off- 
set, which are the offsets that result in maximum-ampli- 
tude I and Q signals. These optimal offset values are io 
then stored in the associated signal channel DDC. 
[0026] At step 404 a selected phase increment is 
stored in the calibration DDC 228, and the I and Q data 
values are obtained. At step 406 an average maximum 
peak-to-peak value for I component is calculated. In the 15 
present embodiment 32 samples are used, but, depend- 
ing upon the application, either more or less samples 
may be used. If the I value calculated at step 406 is 
greater to or equal to a previously obtained average 
maximum peak-to-peak value, then a MAX I 
(f^)peak-to.peak Variable is set equal to the current I value, 
and the corresponding phase value is also saved at step 
410. Otherwise, at step 41 2, the MAX l(n) variable is set 
equal to the previously obtained average maximum 
peak-to-peak value of I, and the corresponding phase 25 
value is saved. 

[0027] This procedure is repeated at step 402 until all 
I phase increments have been tested. The resulting 
MAX l(n) value is then stored into the corresponding 
DDC chip associated with the current channel. The pro- 30 
cedure for determining the Q phase shift value, shown 
in Figure 5, is exactly the same as the procedure for de- 
termining the I phase shift value, as shown in Figure 4. 
These algorithms which provide the calculated phase 
offsets are utilized to move the starting point of the DDC 3S 
208 oscillator signal until it is aligned with the received 
carrier signal, as it enters the DDC 208, to keep the sig- 
nals in phase. 

[0028] Figure 8 is a flowchart illustrating the l/Q bal- 
ance algorithm (step 31 6 of Figure 1 ). This procedure is 40 
necessary in order to normalize the maximum I and Q 
vectors and thus to minimize the errors in phase angle 
calculations. At step 700, the current I and Q gain set- 
tings for the current signal channel are read and the I 
and Q data samples are obtained for the current chan- 45 
nel. At step 704. the I and Q values are compared. If I 
is greater than Q, then a new I gain is calculated at step 
706. The new I gain value is made equal to the current 
I gain value multiplied by the quotient [Q/l]. If, however, 
I is not greater than Q, then a new Q gain value is cal- so 
culated at step 708. The new Q gain value is made equal 
to the current Q gain value multiplied by the quotient [1/ 
O]. The determined gain values are then re-written to 
the DDC 208 of the current signal channel. 
[0029] In the preferred embodiment, the control signal 55 
bus and the data bus for the DDCs 208 are separate. 
The control signal path between the DDCs 208 and the 
receiver card controller 226 is a separate bi-directional 



parallel bus. The data signals are output onto a separate 
serial bus. By using different buses to move data into 
and out of the DDCs 208, the background calibration 
operations are processed in parallel to the high speed 
signal data process flow. As a result, both high speed 
signal processing and background calibration process- 
ing can occur at the same time. In another embodiment, 
both signal and control information share the same par- 
allel bus, and, due to throughput limitations on the bus, 
the calibration routines are run only during startup or 
when there is a break in nomfial system operations. 
[0030] IHowovor the described calibration technique 
only performs the fine tuning of the receiver card 122 
components because its benefits are dependent on the 
parameters provided to the laser modules 116 during 
the initialization time. Thus, there is a need for an im- 
proved calibration method usable with an inboard re- 
ceiver with digital demodulation circuit and fiber optic 
interferometric sensor arrays. Specifically, there is a 
need for an algorithm for coarse calibration of laser sys- 
tems, which can be performed during the system initial- 
ization time. 

[0031] It is thus an object of the present invention to 
provide an apparatus and method for comprehensive 
calibration of laser modules in a signal processing sys- 
tem. 

[0032] This and other objects of the present invention 
may be met by a coarse calibration technique for laser 
modules of a data acquisition module in a signal 
processing system. This coarse calibration technque is 
based on phase modulation of each laser module output 
signal, controlled by a CPU. 

[0033] An optimal value of a phase modulation coef- 
ficient is calculated for each laser and used as the am- 
plitude of a cosine wave signal for the laser module car- 
rier frequency for phase modulation of the laser output 
signal. The phase modulation coefficient is determined 
from a global average l/Q ratio, calculated for alt chan- 
nels of the selected laser module, as an average of all 
l/Q ratios of in-phase and quadrature phase angle com- 
ponents of the signals returned from the sensors. The 
phase modulation coefficient for the selected laser mod- 
ule is determined empirically or by dividing the previous- 
ly determined phase modulation coefficient by the pre- 
viously determined gbbal average l/Q ratio. 
[0034] The exact nature of this Invention, as well as 
its objects and advantages, will become readily appar- 
ent from consideration of the following specification as 
illustrated in the accompanying drawings, in which like 
reference numerals designate like parts throughout the 
figures thereof, and wherein: 

Figure 1 is a high-level block diagram of a signal 
processing system with digital demodulatk)n circuit 
used with the present invention; 
Figure 2 is a block diagram illustrating a preferred 
configuration for an inboard receiver with laser 
nrK)dules driving sensor arrays; 
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Figure 3 is a detailed hardware block diagram ot the 
receiver of Figure 1 ; 

Figure 4 is a flowchart Illustrating the operation of 
the system controller; 

Figure 5 Is flowchart illustrating the phase shift al- 
gorithm for the In-phase signal (I data); 
Figure 6 is a flowchart illustrating the phase shift al- 
gorithm for the quadrature phase signal (Q data); 
Figure 7 Is a flowchart illustrating the functional flow 
of the polarization diversity detector (PDD) selec- 
tion procedure: 

Figure 8 Is a flowchart Illustrating the l/Q balance 

algorithm; 

Figure 9 is a block diagram illustrating a preferred 
configuration for obtaining coarse calibration of a 
single laser module, according to the present inven- 
tion: 

Figure 10 is a flowchart illustrating the steps of the 
coarse calibration algorithm, according to the 
present invention; and 

Figure 11 is a flowchart illustrating the steps of the 
phase modulation coefficient determination module 
of the coarse calibration algorithm, according to the 
present invent bn. 

[0035] The following descriptbn Is provided to enable 
any person skilled in the art to make and use the inven- 
tion and sets forth the best modes contemplated by the 
inventor for carrying out the invention. \^rbus modifi- 
cations, however, will remain readily apparent to those 
skilled in the art, since the basic principles of the present 
invention have been defined herein specifically to pro- 
vide a method and apparatus for coarse calibration of 
laser modules in a data acquisition and signal process- 
ing system, based on phase modulation. 
[0036] The coarse calibration algorithm of the present 
Invention is performed on the laser modules 116 of the 
inboard receiver 102 system with digital demodulation 
circuit and fiber optic Interferometric sensor arrays 104, 
as shown in Figure 1 . Figure 2 and Figure 3. Preferably, 
the coarse calibration of the laser modules 116 is con- 
trolled by the CPU 1 30 and is performed only at the sys- 
tem initialization time, in cooperation with the fine tuning 
calibration method. 

[0037] At the initialization time, some software code 
of the coarse calibration algorithm of the present inven- 
tion, including the data input/output protocol routines, is 
downloaded from the flash memory 1 40 of the CPU 1 30 
to other cards of the inboard receiver 102, such as the 
receiver card 122 and the frequency synthesizer 142. 
The Frequency Synthesizer 142 Includes a processor 
with the Laser Module Supervisor which sets the laser 
modules 116 to proper operating optical frequencies 
and laser module carrier frequencies. The CPU 1 30 has 
the rest of the software code with the coarse calibration 
algorithm. The CPU 1 30 is interfaced with a system con- 
trol panel which is operated by a human operator One 
of the commands entered by the operator is the inboard 



receiver commence initialization command. The opera- 
tor receives the system status report from the CPU 1 30. 
[0038] The download of the calibration code to the fre- 
quency synthesizer 1 42 is accomplished via an RS-232 

5 serial line 1 41 . The CPU 1 30 Viy4E backplane Interface 
139 is attached via a VfviE bus 1 38 to the receiver card 
VME interface 143, using an IDMA port on the receiver 
card controller ADSP-2181 card. The VME bus 138 Is 
used for download from the CPU 130 to the receiver 

10 card 122 over the VME backplane interface 139. 

[0039] The initialization of the lasermodules 116 isac- 
complishod when the Turn on Laser command is sent 
from the CPU 1 30 to the frequency synthesizer 1 42. For 
each laser module 116. the set of laser coefficients Is 

IS read from the laser module 116 EEPROM and sent to 
the frequency synthesizer 142. The laser coefficients 
contain data about the allowed and preferred frequency 
and temperature range for each laser 117. After the ver- 
ification of the laser modules 116 status, which should 

20 have proper voltages and currents, is received in the 
CPU 130 from the frequency synthesizer 142. the star- 
tup of the internal CPU 130 timer Is used for the warm- 
up of the lasers 117. 

[0040] Each laser 117 powers 56 sensor channels of 
2S the array panel 104. The laser 117 crystal is heated to 
the designated temperature level, provided by a signal 
from the frequency synthesizer 1 42, In order to produce 
a laser optical output sine wave signal with the required 
"color'. The carrier frequency Wq also has to be provided 
30 by the frequency synthesizer 142 to the corresponding 
phase modulator 11 9a - 1 1 9h, in the form of a modulat- 
ing cosine wave signal 99. to modulate the phase of the 
laser optical output signal 118. 
[0041] Figure 9 is a block diagram illustrating a pre- 
ss ferred configuration for the coarse calibration of a single 
laser module 116, according to the present Invention. 
The frequency synthesizer 142 controls the operation of 
each laser 1 1 7 by providing, on the amp/f req control line 
146 of Figure 1 , the input cosine wave control signal 99 
40 with the laser frequency Wq (cos Wot) to the correspond- 
ing phase modulator 1 1 9, to modulate the phase of the 
laser optical output signal 118'. In addition, the frequen- 
cy synthesizer 142 provides each laser controller card 
98, and thus the corresponding phase modulator 119. 
45 with the appropriate phase modulation amplitude value 
"C, to control the phase of each laser optical output sig- 
nal 118 in the coarse calibration loop of the Inboard re- 
ceiver 102. 

[0042] The phase modulation amplitude value "C" is 
so calculated in the CPU 130, according to the algorithm 
described below and shown in Figures 11 and 12. 
Therefore, the signal provided to the phase modulator 
119 (C*cos WqI) is used to modulate the phase of the 
corresponding laser output signal 118 to provide the out- 
55 put light signal 1 1 8 with the intensity 1 = A + B'*cos(C*cos 
Wot). The phase of this signal is further phase modulated 
in each interferometric sensor 95 of the array panel 1 04, 
as described in the above-mentioned article by Dan- 
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dridge et al., to obtain the return signal 120 with the in- of tolerance are discarded. At step 1006 it is verified 
tensity I = A + B^cos(C^cos WqI + ^ (t)). Each Interfero- whether the new global average I/O ratio is better than 
metric sensor 95 introduces the change in intensity fronn the previous global average l/Q ratio. If it is determined 
B' to 8, dependent on the mixing efficiency of the inter- that it is better and that the obtained value is close 
ferometric sensor 95, and the phase angle (t) = wt + v s enough to 1 .0 (within appropriate tolerance level pres- 
et), entiy set to O.Oi ), the laser status is set to the calibrated 
[0043] The I and Q components of the phase angle ♦ flag at step 1008 and the routine exits. If it is not better, 
(t) and their ratio l/Q are determined and optimized for a new optimum value of the phase modulation coeffi- 
each sensor channel at the appropriate receiver card cient for the selected laser is determined at step 1010 
1 22 and receiver card controller 226, as described in the io and it is sent to the selected laser at step 1 01 2. 
prior art. They are used both in the fine tuning algorithm [0047] At step 1010. the new optimum value of the 
of the prior art, performed on the recoivGr card 1 22, and phase modulation coefficient for the selected laser may 
in the coarse tuning algorithm of the present invention be determined empirically or by calculation. In the pre- 
performed in the whole inboard receiver 102. ferred embodiment of the present invention this value 
[0044] Figure 10 is a flowchart illustrating the steps of is may be obtained by an appropriate mathematical algo- 
the coarse calibration algorithm of the present invention, rithm performed in the CPU 1 30. The new value of the 
performed under control of the CPU 1 30 at the system phase modulation coefficient, lor example, may be cal- 
initialization time. At step 902 the next uncalibrated laser culated by dividing the old value of the phase modulation 
is selected. After H is determined in step 904 that the coefficient by the old value of the global average l/Q ra* 
selected laser is properly started, at step 906 it is verified so tio. 

whether this is the first calibration attempt for the select- [0048] Therefore, according to the steps of the coarse 
ed laser 117. If it is so, the stored optimum phase mod- calibration algorithm of the present Invention the CPU 
ulation coefficient value is sent to the selected laser, via 1 30 determines the optimum value of the phase modu- 
the frequency synthesizer 1 42. at step 908. Othen^/ise. lation coefficient for each laser module 116. This coeffi- 
the phase modulation coefficient determinatbn routine 2S cient is actually the optimum value of the phase modu- 
of Figure 10 is invoked at step 910. The steps of the lation amplitude C of the cosine wave signal 99 supplied 
coarse calibration algorithm are repeated iteratively for to the phase modulator 11 9 by the frequency synthesiz- 
each laser, as long as is needed to obtain the l/Q ratio er 142, as shown in Figure 8. It is possible that each 
for each channel as close as possible to 1.0. at which sensor in the array panel 1 04 receives a slightly different 
time it is proclaimed at step 91 2 that ail lasers controlled 30 signal from the assigned laser 1 1 7, so it is impossible to 
by that CPU are calibrated. At step 91 4 the optimum val- provide a perfect value for the phase modulation ampli- 
ues of the phase modulation coefficients for all laser tude C which would be ideal for all sensors, but only a 
modules controlled by that CPU are stored in the CPU value optimal for most sensors. Therefore, each receiv- 
1 30 for the next initialization time. er card 1 22 channel has to be fine tuned, as described 
[0045] Figure 1 1 is a flowchart illustrating the steps of 3S in the above mentioned prior art patent applications, 
the phase modulation coefficient determination module [0049] In one embodiment of the present invention 
ofthecoarsecalibratlon algorithm, performedforthese- each CPU 130 has a memory map from which It is de- 
lected laser according to the present invention. At step termined which receiver cards 122, controlled by that 
1002 a command is sent to the appropriate receiver card CPU, are associated with the selected laser module 
(s) 1 22 receiver card controller 226 to collect the l/Q ratio 40 ^^s. Each CPU 1 30 preferably controls eight lasers. Ide- 
value for each channel assigned to the selected laser, ally, each CPU 1 30 assists all 56 channels of the laser 
calculated from the average maximum Ipeak-io-peak module 116 when they are all received on the same re- 
average maximum Qpeak-to-peak values, and obtained ceiver card 122. However, for practical reasons, in the 
according to the algorithm of Figure 3. The l/Q ratio is Implementation of the present invention each receiver 
calculated for each channel as a ratio of an average 45 card 122assistsonly 28 channels of each laser module 
maximum lvalue (i.e., when Q is near 0) divided by the 116. Therefore, each array panel 1 04 sensor channels 
average maximum Q value (i.e. , when I is near 0). where are divided between two CPUs and the memory map in 
the average I and Q values are taken from a predeter- this embodiment defines which CPU coarse calibration 
mined number of samples for that channel. Therefore, module would be used to perform the calibration of the 
the l/Q values calculated during the fine tuning calibra- so whole array panel 1 04, based on data from the 28 chan- 
tion at the system initialization time, at the appropriate nels assigned to that CPU. It is conceivable, however, 
receiver card 1 22, are sent to the CPU 1 30 via its VME that in another embodiment each CPU would calculate 
backplane interface 1 39. its own global average l/Q ratio and that the value closer 
[0046] At step 1004, performed in the CPU 130, the tol.O.ortheaverageofthetwovalues, would be chosen 
collected values of l/Q ratio for all channels controlled ss to determine the phase modulation amplitude C for that 
by the selected laser are added, in order to determine laser. 

a global average l/Q ratio for all the channels controlled [0050] Those skilled in the art will appreciate that var- 

by the selected laser. The sensor values which are out bus adaptations and modifications of the just-described 
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preferred embodirr^ents can be configured without de- 
parting from the scope and spirit of the invention. There- 
fore, it is to be understood that within the scope of the 
appended claims, the invention may be practiced other 
than as specifically described herein. For example, in 
addition to the above mentioned undenvaier applica- 
tions the present invention is also applicable to commer- 
cial and military seismic applications. 



Claims 

1. A calibration method for laser modules in a signal 
processing system, comprising the steps of: 

selecting a next laser module to be calibrated; 
determining if the selected laser module has 
been calibrated before; 

if not calibrated before, sending a saved phase 
modulation coefficient to the selected laser 
module to be used for phase modulation of the 
laser module output signal; 
if calibrated before, determining a new phase 
modulation coefficient to be used for phase 
modulation of the laser module output signal; 
repeating the above steps for each laser mod- 
ule in the signal processing system until all la- 
ser modules are calibrated; and 
saving the new phase modulation coefficients 
for all laser modules. 

2. The method of claim 1 , wherein the step of deter- 
mining the new phase modulation coefficient further 
comprises the steps of: 

obtaining ratios of in-phase and quadrature 
phase angle components (l/Q) ratios for alt 
channels of the selected laser module; 
calculating a new global average l/Q ratio for 
all channels of the selected laser module as an 
average of all said l/Q ratios; 
comparing the new global average l/Q ratio with 
the previously determined global average l/Q 
ratio; and 

if the new global average l/Q ratio is closer to 
1 than the previously determined global aver- 
age l/Q ratio and within a predetermined toler- 
ance, giving the selected laser module a cali- 
brated status, or 

if the new global average l/Q ratio is not closer 
to 1 than the previously determined global av- 
erage l/Q ratio or not within the predetermined 
tolerance, determining a new phase modula- 
tion coefficient for the selected laser module; 
and 

sending the new phase modulation coefficient 
to the selected laser module to be used for 
phase modulation of the laser rrKxjule output 
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signal as the amplitude of a cosine wave signal 
with the laser nrxxiule carrier frequency. 

3. The method of claim 1 or claim 2, wherein the step 
of determining the new phase modulation coeffi- 
cient for the selected laser module performs deter- 
mination by an empirical method. 

4. The method of any of claims 1 to 3, wherein the step 
of determining the new phase modulation coeffi- 
cient for the selected laser module further compris- 
es the step of: 

obtaining the new phase modulation coeffi- 
cient by dividing the previously determined phase 
modulation coefficient by the previously determined 
gbbal average l/Q ratio. 

5. The method of any of claims 1 to 4, wherein the step 
of obtaining ratios of in-phase and quadrature 
phase angle components (l/Q ratios) for all chan- 
nels of the selected laser nruxJule further comprises 
the steps of: 

calculating an average maximum value for an 
in-phase component of a phase angle received 
for the channel; 

calculating an average maximum value for a 
quadrature component of the phase angle re- 
ceived for the channel; 

dividing the calculated average maximum val- 
ue for the in-phase component of the phase an- 
gle with the calculated average maximum value 
for the quadrature component of the phase an- 
gle; and 

repeating the above steps for each channel of 
the selected laser module. 

6. The method of any of claims 1 to 5, wherein the step 
of calculating the average maximum value for the 
in-phase component of the phase angle further 
comprises the steps of: 

obtaining a plurality of maximum peak-to-peak 
values of the in-phase component, taken when 
the quadrature component is near or at zero; 
and 

calculating the average maximum value for 
said in-phase component from the plurality of 
maximum peak-to-peak values of the in-phase 
component. 

7. The method of any of claims 1 to 6, wherein the step 
of calculating the average maximum value for the 
quadrature component of the phase angle further 
comprises the steps of: 

obtaining a plurality of maximum peak-to-peak 
values of the quadrature component, taken 
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when the in-phase component Is near or at ze- 
ro; and 

calculating the average maximum value lor 
said quadrature component trom the plurality of 
maximum peak-to-peak values of the quadra- s 
ture component. 

8. An apparatus for calibrating a signal processing 
system, the apparatus comprising: 

10 

at least one laser module sending an output 
light signal to an associated array of sonsor 

channels; 

at least one receiver for receiving return signals 
from each of the sensor channels, the return is 
signals having an in-phase component and a 
quadrature component of a phase angle; 
a system processor for receiving the return sig- 
nals from the sensor channels and determining 
a phase modulation coefficient for the laser 20 
module; and 

a phase modulator tor receiving a cosine wave 
signal at the carrier frequency and the ampli- 
tude equal to the phase modulation coefficient, 
the cosine wave signal used for phase modu- 2S 
lation of the laser module optical output signal. 

9. The apparatus of claim 8 wherein said phase mod- 
ulation coefficient is dependent on a global average 

l/Q ratio, calculated for all the sensor channels of 30 
the laser module, wherein said global average l/Q 
ratio is calculated as an average of the ratios of the 
In-phase and quadrature phase angle component 
(l/Q) received for all sensor channels of the laser 
module. ^ 

10. The apparatus of claim 8 or claim 9, wherein said 
phase modulation coefficient is determined empiri- 
cally. 

40 

11. The apparatus of claim 8 or claim 9, wherein said 
phase modulation coefficient is calculated by divid- 
ing the previously determined phase modulation co- 
efficient by the prevbusly determined global aver- 
age l/Q ratio. ^ 

12. The apparatus of any of claims 8 to 11, wherein 
each said ratio of in-phase and quadrature phase 
angle component (l/Q ratio) is determined in the ap- 
propriate receiver, which cateulates an average so 
maximum value for an in-phase component of a 
phase angle and divides it with an average maxi- 
mum value for a quadrature component of the 
phase angle. 

55 

13. The apparatus of claim 12, wherein the average 
maximum value for the in-phase component of the 
phase angle is calculated as an average of a plu- 



rality of maximum peak-to-peak values of the in- 
phase component, taken when the quadrature com- 
ponent is near or at zero, and 

the average maximum value for the quadra- 
ture component of the phase angle is calculated as 
an average of a plurality of maximum peak-to-peak 
values of the quadrature component, taken when 
the in-phase component is near or at zero. 
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